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We report an observation of the decay — > (fxf>K ± and evidence for B° — > <fxf>K°. These results 
are based on a 414 fb" 1 data sample collected with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider operating at the T(4S) resonance. The branching fractions for these decay 
modes are measured to be B(B ± -> = (3.2±JJj ± 0.3) x 10~ 6 and B{B° -> HK°) = 

(2.3±J;? ± 0.2) x 10" 6 for <j>(f> invariant mass below 2.85 GeV/c 2 . The corresponding partial rate 
asymmetry for the charged B mode is measured to be Acp(B ± — > 4>4>K ± ) = O.OIIq' ]_t 

±0.02. We 

also study the decays B — > J/ipK ± and B — > rjcK^ , where the J/^ and ?7 C decay to final states 
with four charged kaons. We find Acv(B ± (p(f>K ± ) with the </></> candidates within the rj c mass 
region is 0.15lto'i7 i 0.02, consistent with no asymmetry. 



PACS numbers: 13.25.Ft, 13.25.Hw, 14.40.Nd 

Evidence of charmless B — > (f><j>K decays has been 
reported by the Belle collaboration using 85 x 10 6 BB 
pairs In the Standard Model (SM), this decay chan- 
nel requires the creation of an additional final ss quark 
pair in a b — > sss process, such as B — > 0if . Therefore, 
the study of P> — + </xf>K provides useful information for 
understanding quark fragmentation in B decays. More- 
over, our previous study also reported results for the de- 
cays B — > J /ij)K(r) c K) with the J/ip(r] e ) decaying into 
four kaons in the final state, which can proceed with <f> 
mesons in the intermediate state. It has been suggested 
that large direct CP violation up to 40% is possible in 
B — ► rj c K — ► cf><j>K if there is a sizable b — > s transition 
with a non-SM CP-violating phase that interferes with 
the decay amplitude via the rj c resonance [2|. 

Recently the BaBar collaboration has reported results 
of a study of B -> ^K 0. The branching fraction 
for i? ± — > <f>qbK^ that they obtained is around three 
times larger than our previous measurement Here 
we present improved measurements of B — > (/></>_ftT decays 
with not only larger statistics but also proper consid- 
eration of the non-resonant K + K~ contribution under- 
neath the 4> resonance. The analysis is based on a data 
sample of 414 ftT 1 containing 449 xl0 6 BB pairs. The 
data were collected with the Belle detector at the KEKB 
asymmetric-energy e + e~ (3.5 on 8 GeV) collider [H op- 
erating at the T(45) resonance. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters ( ACC) , a barrel-like ar- 
rangement of timc-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect 



K\ mesons and to identify muons (KLM). The detector 
is described in detail elsewhere [5J. Two inner detector 
configurations were used. A 2.0 cm radius beampipe and 
a 3-layer silicon vertex detector (SVD1) were used for 
the first sample of 152 xl0 6 £?-B pairs, while a 1.5 cm 
radius beampipe, a 4-layer silicon detector (SVD2) and 
a small-cell inner drift chamber were used to record the 
remaining 297 xlO 6 BB pairs 0. 

Charged kaons are required to have impact parame- 
ters within ±2 cm of the interaction point (IP) along the 
z-axis (antiparallel to the positron direction) and within 
0.2 cm in the transverse plane. Each track is identified 
as a kaon or a pion according to a Kjix likelihood ratio, 
1Z(K/tt) — Ck/{Ck + At), where Ck/£-k is the likeli- 
hood of kaons/pions derived from the responses of TOF 
and ACC systems and the energy loss measurements from 
the CDC. The likelihood ratio is required to exceed 0.6 
for kaon candidates; within the momentum range of in- 
terest, this requirement is 88% efficient for kaons and has 
a misidentification rate for pions of 8.5%. Neutral kaons 
are reconstructed via the decay Kg — > tt + tt~ and have an 
invariant mass 0.482 GeV/c 2 < M w+V - < 0.514 GeV/c 2 
(±4cr mass resolution). The 7r + 7r _ vertex is required to 
be displaced from the IP and the flight direction must be 
consistent with a K s that originated from the IP. The 
required displacement increases with the momentum of 
the Kg candidate. 

B meson candidates are reconstructed in the five-kaon 
final state. Two kinematic variables are used to distin- 
guish signal candidates from backgrounds: the beam- 
energy constrained mass M hc = \J El eara - |P reC on| 2 and 
the energy difference AE = E recoIi - Pbcam, where E heam 
is the beam energy, and P rcC on and P r ccon are the re- 
constructed energy and momentum of the signal candi- 
date in the T(4S') rest frame. The resolution of M^ c is 
approximately 2.8 MeV/c 2 , dominated by the beam en- 
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ergy spread, while the AE resolution is around 10 MeV. 
Candidates with five kaons within the region \AE\ < 
0.2 GcV and 5.2 GeV/c 2 < M bc are selected for fur- 
ther consideration. The signal region is defined as 5.27 
GeV/c 2 < M bc < 5.29 GeV/c 2 and \AE\ < 0.05 GeV. 

The dominant backgrounds are e + e~ — > qq (q = 
U, d, s, c) continuum events. Event topology and B fla- 
vor tagging are used to distinguish the jet-like continuum 
events and the spherically distributed BB events. Seven 
event-shape variables are combined into a single Fisher 
discriminant Q- The Fisher variables include the an- 
gle between the thrust axis of the B candidate and the 
thrust axis of the rest of the event (cos 9t), five modified 
Fox- Wolfram moments 0] , and a measure of the momen- 
tum transverse to the event thrust axis (S±) 9]. Two 
other variables that are uncorrelated with the Fisher dis- 
criminant and help to distinguish signal from the con- 
tinuum are cos 9 b , where 9 b is the angle between the 
B flight direction and the beam direction in the e + e~ 
center-of-mass frame, and Az, the z vertex difference be- 
tween the signal B candidate and its accompanying B. 
We form signal and background probability density func- 
tions (PDFs) for the Fisher discriminant, cos 9b and Az 
using the signal Monte Carlo (MC) events and sideband 
data (5.2 GeV/c 2 < M hc < 5.26 GeV/c 2 ), respectively. 
The products of the PDFs for these variables give signal 
and background likelihoods £$ and Cbg f° r each candi- 
date, allowing a selection to be applied to the likelihood 
ratio K = C s / {Cs + C B g)- 

Additional background discrimination is provided by 
the quality of the B flavor tagging of the accompany- 
ing B meson. The standard Belle flavor tagging pack- 
age gives two outputs: a discrete variable indicating 
the flavor of the tagging B and dilution factor r, which 
ranges from zero for no flavor information to unity for 
unambiguous flavor assignment. The continuum back- 
ground is reduced by applying a selection requirement 
on the variable 1Z for events in each r region according 
to the figure of merit defined as Ns/^/Ns + Nbg, where 
Ns denotes the expected <fxf>K signal yield based on MC 
simulation and the branching fraction reported in our 
previous measurements, and Nbg denotes the expected 
qq yields from sideband data. This requirement removes 
(61-81)% of the continuum background while retaining 
(80-92)% of the signal, and depends on the decay chan- 
nel (cjxfiK or (f>cj)K Q ) and the SVD configuration during 
the measurement (SVD1 or SVD2). Backgrounds from 
other B decays are investigated using a large MC sample 
and are found to be negligible after the 1Z requirement. 

The signal yields are extracted by applying an un- 
binned extended maximum likelihood (ML) fit to the 
events with M bc > 5.2 GeV/c 2 and \AE\ < 0.2 GeV. 
For the (jx^K^ mode, we simultaneously obtain the yield 
and the partial rate asymmetry Acp using the likelihood, 
defined as: 
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c = e - {Ns+NBG) JJ(£ ±[1 - q t . A^N.Pl), (1) 
» 3 

where i is the identifier of the i-th event, j indicates sig- 
nal (S) or background (BG), P is the two-dimensional 
PDF of Mbc and AE, and q indicates the B meson fla- 
vor, + 1 for B + and —1 for _B~, respectively. For neu- 
tral B events, the factor |[1 — ■ Acp] in Eq. ((T|) is 
replaced by 1. The Mb c PDFs are modeled by a Gaus- 
sian function for signals and an ARGUS function [Tljj for 
the continuum, while a Gaussian is used to describe the 
signal AE and a second-order Chebyshev polynomial is 
used for the background AE distribution. The param- 
eters of the PDFs are determined using high-statistics 
MC samples and sideband data for signal and back- 
ground shapes, respectively. The signal PDFs are cal- 
ibrated by comparing the Mb c and AE distributions of 
the B+ ->■ D°{K+n-TT-'K + )TT + data with the MC expec- 
tation. 

We search for charmless B — > (fxpK decays by requir- 
ing the cfxj) invariant mass (M^) to be less than 2.85 
GeV/c 2 , the region below charmonium threshold. Can- 
didate <t> mesons are identified by requiring the invariant 
masses of K + K~ pairs (M K + K -) to be in the range 1.0 
GeV/c 2 to 1.04 GeV/c 2 (±4.6<r). Figure CD shows the 
Afbc and AE projections with the fit curves superim- 
posed. Clear signals appear in both B^ and B° modes 
with signal yields of 37.0rJ; an< ^ t-&-2.5> respectively. 
Although K + K~ candidates are required to lie in the 
(j) mass region, non-0 backgrounds may also contribute. 
Figure^a) shows the M^ +if _ vs. M^ c+K _ distributions 
for (K + K~K + K~)K ± candidates in the signal region, 
where the two K + K~ pairs are required to have invari- 
ant masses less than 1.2 GeV/c 2 . Events in the two <j> 
bands are used to estimate the B ± — > <j)K + K~ con- 
tribution. Figure [2jb) shows the B signal yields [13] as a 
function of the K + K~ invariant mass after requiring the 
other K + K~ pair to have a mass in the <fi mass region. 
The B signal yields are fitted with a threshold function 
in the region 0.98 GeV/c 2 < M K+K - < 1.20 GeV/c 2 , 
excluding the (j> mass region (1.00 GeV/c 2 < M K + K ~ < 
1.04 GeV/c 2 ). The size of the non-0 contribution is es- 
timated by interpolating the B yields in the <f> sideband 
region to the <fi mass region. Since events in the two 
4> bands contain both true (j) mesons and non-resonant 
K + K~ pairs, the area underneath the <f> mass region in 
Fig. [2jb) also includes the <\>K + K~K ± contribution and, 
due to combinatorics, a double counted contribution from 
the non-resonant 5K component. The contribution of 
B — > 5K is estimated by extrapolating the B signal yield 
in the upper right corner of the dashed region in Fig. [2](a) 
to the 4> mass region. The fraction of non-<j)4>K events in 
the 4> mass region as obtained from both contributions 
is thus (7 ± 4)%. The same procedure is applied to the 
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(f>(f>K sample; here we obtain a fraction of (7 ± 9)%. 
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FIG. 2: (a) The distribution of M 1 K+K _ vs 
the K+K-K+K'K* candidates in the M hc - AE signal 
box with M K+K - < 1.2 GeV/c 2 . The two K + K~ bands 
indicate the mass region (1.0 GeV/c 2 < M K+K - < 1.04 
GeV/c 2 ). The horizontal and vertical dashed lines are lo- 
cated at M K + K - — 1.05 GeV/c 2 . The rectangle on the up- 
per right is the 4>4> sideband region; events in this region are 
used to estimate the non-resonant B — > 5K contribution, (b) 
B signal yield as a function of the M K + K - after requiring 
the other K + K~ pair to have a mass in the <j> mass region. 
A threshold function is used to fit the data; events with 1.0 
GeV/c 2 < M K+K - < 1.04 GeV/c 2 are excluded when the fit 
is performed. 



FIG. 1: Distributions of M^c and AE for the decay modes 
B ± -» (a,b) and B° -> 4>cj>K° (c,d), with 00 in- 

variant mass less then 2.85 GeV/c 2 . The open histograms 
are the data, the solid curves show the result of the fit, the 
dash-dotted curves represent the signal contributions and the 
dashed curves show the continuum background contributions. 

Table Q] summarizes the (f>(j)K results after subtract- 
ing the non-<fi(f>K contribution. Signal efficiencies are ob- 
tained by generating c/xf)K MC events, where the same 
< 2.85 GeV/c 2 requirement is applied. System- 
atic uncertainties in the fit are obtained by performing 
fits in which the signal peak positions and resolutions 
of the signal PDFs are successively varied by ±1ct. The 
quadratic sum of each deviation from the central value 
of the fit gives the systematic uncertainty of the fit. For 
each systematic check, the statistical significance is taken 

as y/ — 2 ln(£feeddown/£max), where £feeddown and £max 

are the likelihoods at the expected non-iptpK yields and 
the best fit, respectively. The change in significance that 
arises from uncertainties in the signal PDFs is negligi- 
ble (< 1%). We choose the significance calculated after 
increasing the non-00i4T yield by its la statistical un- 
certainty as our significance including systematic uncer- 
tainty The numbers of B + B~ and B a B° pairs are as- 
sumed to be equal. 

The systematic uncertainty resulting from the TZ re- 
quirement is studied by checking the data-MC efficiency 
ratio using the B + — > D°(-^ K + tt~ ir~ tt + )tt + sample. 
The corresponding systematic error is 2.7-2.8% and again 
depends on the decay channel and SVD geometry. The 
systematic errors on the charged track reconstruction are 



estimated to be around 1% per track using partially re- 
constructed D* events. Therefore, the tracking system- 
atic error is 5% (five tracks) for the 00if ± mode and 
4% for the (fi(f>K mode (excluding Kg reconstruction). 
The kaon identification efficiency is studied using samples 
of inclusive D*+ -> D 7r+,D° -> K~-k+ decays. The 
Kg reconstruction is verified by comparing the ratio of 
D+ -> K s tt+ and D+ -> K-tt + tt+ yields. The resulting 
Kg detection systematic error is 4.9%. The uncertainty 
in the number oiBB events is 1.4%. The final systematic 
error is obtained by summing all correlated errors linearly 
and then quadratically summing the uncorrelated errors. 

After subtracting the non-cfxpK contribution, the 
branching fractions for charmless B — > (jxpK decays are 
B(B ± -► (faK*) = (3.2±g:| ± 0.3) x 1CT 6 with a 9.5a 
significance and B(B° -> 00X ) = (2.3±J:? ± 0.2) x 10~ 6 
with a 4.7cr significance. The measured charge asymme- 
try for B ± -> 00/^ decay is 0.01+°,;^ ± 0.02. The first 
error is statistical and the second is systematic. 

It is of interest to search for possible 00 resonances 
above charmonium threshold. Figure [3^a) shows the B 
signal yield divided by the bin size as a function of 
if the < 2.85 GeV/c 2 requirement is not applied. 

There is no enhancement in the high 00 mass region 
except for the r\ c peak near 3 GeV/c 2 . Reference [| 
suggests the possibility of a large CP asymmetry, which 
could arise from the interference between B^ — > 00i4T ± 
and B^ — ► rj c (—> 00)i"C ± decays. Events with 00 invari- 
ant mass within ±40 MeV/c 2 of the nominal r\ c mass are 
selected to investigate this asymmetry. The measured 
CP asymmetry is 0.15lg 17 ± 0.02, which is consistent 
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TABLE I: Mode, yield, efficiency including secondary branching fractions, branching fraction for B — > <p<pK and related 
charmonium decays. 

Mode Yield Efficiency(%) B(10" 6 ) 

#± _^ ^jf± (M 00 < 2.85 GeV/c 2 ) 34.2^;| 2ll 3.2+^ ± 0.3 

B° -> <txt>K° (M^ < 2.85 GeV/c 2 ) 7.3±|;2 0.69 2.3± ;? ± 0.2 

B ± — > r\ c K ± , r\ c — > ^ 27.9±gj 2.72 2.3 ± 0.6 ± 0.2 

-> r\ c K ± , r) c -> <j}K+K~ 60.3±^;^ 4.85 2.8±E^ ± 0.2 

B ± -> rfc-K - *, jj c -> 2(JT+if-) 105.7±^:^ 9.93 2.4tg;{j ± 0.2 

-> J/ipK*, J/tp -> (pK+K~ 26.3tjj;? 4 -67 1.3 ± 0.3 ± 0.2 

B ± J/tp^, J/tp — ► 2(K + K~) 36.0± 7 7 l 9.41 0.85±° -\ 8 7 ± 0.10 



with no asymmetry. 

We study possible charmonium states by measuring 
the B yield with M iK between 2.8 GeV/c 2 and 3.2 
GeV/c 2 . Since r\ c and J/tp mesons may decay to (j)K + K~ 
and 2(K + K~), mass fits are performed with and with- 
out the requirement that one or both K + K~ pairs lie in 
the 4> mass region. As shown in Fig. [3J clear rj c and J/tp 
resonances are visible in the (pK + K~ and AK samples 
while only an r\ c peak appears in the (p<p mode. 

We obtain the signal yields for B ± — > ijcK^ and 
B ± -> J/V'^ ± by performing x fits with asymmetric 
errors to the M^, M^k+k- and M^k invariant mass 
distributions, which are presented in Figs. [3][b, c, d). 
The J/tp signal PDF is modeled with a Gaussian func- 
tion while the r\ c PDF is described by a Breit-Wigner 
function convoluted with a Gaussian resolution function, 
which has the same Gaussian width as the J/tp PDF. 
Since sizable signals are observed in the AK mode, the 
parameters are determined using the AK sample and the 
same signal PDFs are then applied to the cpK + K~ and 
(pep samples. The obtained Gaussian width is measured 
to be 4.5+ 1 ; 3 MeV/c 2 . The resulting signal yields are 
summarized in Table Q] The peak positions obtained for 
the r) c and J/tp are 2.979 ± 0.002 ± 0.001 GeV/c 2 and 
3.094 ± 0.001 GeV/c 2 , respectively, consistent with the 
nominal rj c and J/tp masses. The f] c Breit-Wigner width 
is measured to be 29.8±g 2 5 2 ±0.1 MeV/c 2 , where the cen- 
tral value is consistent with the world average [3] and 
the second error is due to the uncertainty of the Gaussian 
width for the mass resolution. The validity of determin- 
ing B signal yields from a constrained x 2 fit with an 
asymmetric error is verified by toy MC. 

For the <pK + K~ and <p<p modes, the non-0 contribu- 
tion is determined from the B signal yields for events 
with one K + K~ pair in the <p sideband region (1.05 
GeV/c 2 < M K+K - < 1.09 GeV/c 2 ) and the AK and 
cpK + K~ masses in the charmonium resonance region, re- 
spectively. We find 3.0^ 2 '4 events in the rj c — ► (p<p mode, 

1+5.4 



6.4_ 4 ; 5 events in the r] c 



(PK+K- mode, and 3.5^^ 
events in the J/tp — > <pK + K~ mode. For the yields of 
these modes, listed in Table HI the corresponding feed- 
down yields have been subtracted. 

Signal efficiencies are determined using signal MC and 
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3: B signal yield as a function of (a,b) M^, (c) 
^ X+Jf _ and (d) Mak- In (a) we use different bin sizes for 
less than 3 GeV/c 2 and greater than 3 GeV/c 2 . The sub- 
set with from 2.8 GeV/c 2 to 3.2 GeV/c 2 is shown in (b). 
The J/tp signal PDF is modeled with a Gaussian function 
while the rj c PDF is described by a Breit-Wigner function 
convoluted with a Gaussian resolution function, which has 
the same Gaussian width as the J/tp PDF. The solid curves 
in (b,c,d) show the results of the fit and the contributions 
(second order polynomial) not from the J/tp and r\ c . 



their systematic uncertainties are similar to what was 
described in the charmless cpcpK part. Systematic uncer- 
tainties in the fitting are obtained by performing fits in 
which the signal peak positions, the resolutions of the sig- 
nal PDF's and the width of the Gaussian resolution func- 
tion convoluted with the r\ c PDF are successively varied 
by ±er. The quadratic sum of all deviations gives the sys- 
tematic error of the fit. The products of the branching 
fractions for various decays are listed in Table |TJ Since 
the probabilities of r\ c and J / tp decays to AK, cpKK and 
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(f>4> final states are measured with poor accuracy, we can 
use the values from Table U as well as the world average 
branching fractions B(B ± r/ c X ± ) = (9.1 ± 1.3) x 10" 4 
and B(B ± -> J/ipK*) = (1.007 ± 0.035) x 10" 3 [H 
to determine independently the corresponding branching 
fractions for the r\ c and J/tp; the results are shown in Ta- 
ble [TTJ The world average values given above are based 
on datasets, which also include Belle measurements fl4j ; 
however, our estimate shows that effects of possible cor- 
relations are negligible. We also provide the ratios of 
branching fractions as shown in Table IIIII The system- 
atic uncertainties are predominately due to the Gaussian 
width for the mass resolution; other systematic uncer- 
tainties either cancel out in the ratios or too small to be 
considered. 



ton collisions. For both the r\ c and J/i[> decays into four 
charged kaons, the <fiK + K~ contribution dominates. 
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TABLE II: The measured branching fractions of secondary 
charmonium decays. 



Decay mode 



B 



Vc -> <f,K + K- 
T) C -> 2(K + K~) 
J/ip -> 4>K + K~ 
J/rp -> 2(K+K- 



(2.5ig;g ±0.4) x 10 ~ 
(3.0 ±0.6 ±0.5) x 10~ 3 

(2.6lo;5 ± 0.4) x 10~ 3 
(1.2 ±0.3 ±0.2) x 10~ 3 



± 1.0) x 10" 



TABLE III: Ratios of branching fractions for the rjc and J/ip 
decays. 



Mode 




Ratio 


B(,, c ->4>0) 




0.961^2° ± 0-02 
1.17l°;f 6 ± 0.01 

i h-7+0.49+0.06 
LAI -0.46-0.08 


B(rj c ^2(X+ K~ 


)) 
) 


B(rj c ^2(i?+ K~ 
B(J/ip^tf>K+ K 


)) 
-) 


B(J/ip^2(K+K- 


")) 



In summary, we have observed the charmless decay 
B ± — ► </></> if ± and evidence for B° — > cfxfiK . We also re- 
port the CP asymmetry of the charged decay and mea- 
surements of other closely related charmonium decays. 
The results are consistent with our previous measure- 
ments [l| and supersede them, and have considerably 
better precision due to the increase in statistics. The 
obtained branching fraction of charmless B ± — > (fxpK^ 
decay is 3.1<r smaller than the measurement by the BaBar 
collaboration For the charmonium decays, we find 
that the decay rj c — > (jxp contributes around 42% and 
24% of the r] c -► cj>K + K~ and r] c -► 2(K+K~) events, 
respectively. The latter is consistent with an early Belle 
measurement using r\ c events produced in two pho- 
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